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Abstract 

In this paper, we consider the higher order Boussinesq (HBq) equation which 
models the bi-directional propagation of longitudinal waves in various con¬ 
tinuous media. The equation contains the higher order effects of frequency 
dispersion. The present study is devoted to the numerical investigation of 
the HBq equation. For this aim a numerical scheme combining the Fourier 
pseudo-spectral method in space and a Runge Kutta method in time is 
constructed. The convergence of semi-discrete scheme is proved in an ap¬ 
propriate Sobolev space. To investigate the higher order dispersive effects 
and nonlinear effects on the solutions of HBq equation, propagation of single 
solitary wave, head-on collision of solitary waves and blow-up solutions are 
considered. 

Keywords: The higher-order Boussinesq equation, Fourier pseudo-spectral 
method. Solitary waves. Head-on collision of solitary waves. Blow-up 
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1. Introduction 

The hrst attempt to describe the solitary waves of Scott Russell’s exper¬ 
iment in 1834, theoretically, was made by Boussinesq 0,3 3 - To model 
the motion of the long waves in shallow water he proposed the Boussinesq 
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( 1 . 1 ) 


We refer the reader to [sf for a more detailed information on the origins of 
Boussinesq equation. As the equation (II.ip turned out to be incorrect in the 
sense of Hadamard, several re-derivations are made. The most common well- 
posed re-derivation is proposed for ion-sound waves and for quasi-continuum 
approximation of dense lattices [^,[^ and is called the ’’improved” Boussinesq 
equation (IBq) 

Utt ~ '^xx "^xxtt (/*(^)) XX' ( 1 . 2 ) 


Unlike the Boussinesq equation (II.ip . the IBq equation is not fully integrable. 
The Boussinesq equation and its re-derivations include the lowest order effects 
of frequency dispersion. To include the higher order effects of dispersion, 
Rosenau proposed the higher order Boussinesq (HBq) equation 


'^XX T Vl^xxtt V‘2'^xxxxtt T (/(^))xx (1"3) 


where rji and 7^2 are real positive constants. In Rosenau’s own words ’’low¬ 
est order dispersion shows some anomalous which can be resolved only if 
higher order discrete effects are included.” A more recent derivation of the 
HBq equation is given by Duruk et ah in to model the bi-directional 
propagation of longitudinal waves in an inhnite, nonlocally elastic medium. 
Eliminating the effects of the higher order dispersion, the HBq equation re¬ 
duces to the IBq equation. 

The local and global well-posedness of the Cauchy problem for the HBq 
equation together with the initial conditions 


u{x,0) = (f){x), Ut{x,0) = tlf^x) (1.4) 

was proved in in the Sobolev space with any s > 1/2. The same 
authors also studied the qualitative properties of a general class of nonlocal 
equations which includes the HBq equation as a special case in [lol |. 

The present study addresses the answers of the following very natural 
questions: How the higher-order dispersive term and the various power type 
nonlinear terms, i.e. f{u) = p > 1, affect the solutions? Do so¬ 

lutions of the HBq equation converge to solutions of the generalized IBq 
equation as r ]2 O’*"? To answer these questions we need an efficient nu¬ 

merical method. We therefore propose a numerical method combining a 
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Fourier pseudo-spectral method for the space discretization and a fourth- 
order Runge-Kutta scheme for time discretization. To the best of our knowl¬ 
edge, there is no numerical study for HBq equation although there are lots 
of numerical studies to solve the generalized IBq equation (see 0 and refer¬ 
ences therein). In Section 2, we review some properties of the HBq equation 
and then we derive the solitary wave solutions. In Section 3, the semi-discrete 
scheme is introduced for the HBq equation and the convergence of the semi¬ 
discrete scheme is proved in an appropriate space. In Section 4, we propose 
the fully-discrete Fourier pseudo-spectral scheme and show how to formu¬ 
late it for the HBq equation. In Section 5, the numerical scheme is tested 
for accuracy and convergence rate. The effect of extra dispersion and the 
nonlinearity on the solutions of the HBq equation by considering various 
problems, like propagation of single solitary wave, head-on collision of two 
solitary waves and the blow-up solutions are discussed in Section 6. 


2. Properties of the higher-order Boussinesq equation 


In this section, we will review the conservation laws of the HBq equation 
and we will derive the solitary wave solution of the HBq equation by using 
ansatz method. 

(a) Conservation laws: Three conserved integrals in terms of U where 
u = Ux for the HBq equation are given as 


hit) 

hit) 

hit) 



dx, 


Ux\Uf T^iUxxt T V^Uxxxxt) dx 


-OO 

f*00 


( 2 . 1 ) 

( 2 . 2 ) 


l(u,f + 2F(U,) + (Urf + dx (2.3) 


with /(s) = F'is). The derivation of the integrals can be found in js], 10 
(b) Solitary wave solution: We use the ansatz method which is the most 
effective direct method to construct the solitary wave solutions of the non¬ 
linear evolution equations. We look for the solutions of the form u = u^) 
where ^ = x — ct — xq under asymptotic boundary conditions. Substituting 
the solution u = m(^) into equation fll.3p and then integrating twice with 
respect to we have 


l)u 


9 9 d'^u 


(2.4) 
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We now look for the solution of the form 


m(^) = Asech'^(i?^). 


(2.5) 


If we substitute the above ansatz into the equation fl2.4p . the solitary wave 
solution of the HBq equation is given in the form 


u{x, t) = A {sech"^ {B{x — cf — Xq)) } , (2.6) 


r]lc^{p -M) {p + 3) (3p 1) 

P-1 

, B = 

1 

1 

to 

_1 

2p2 (p2 + 2p + 5)2 


4^72 (p^ + 2.p + b) 

to 

1 

Apl {p+1) ^ 

_P2 + 2.p + 5)2 _ 

y 


where A is the amplitude and B is the inverse width of the solitary wave. 
Here c represents the velocity of the solitary wave centered at Xq with (? > 1. 


3. The semi-discrete scheme 


3.1. Notations and Preliminaries 

Throughout this section, C denotes a generic constant. We use (•, •) and 
II ■ II to denote the inner product and the norm of dehned by 

{u,v)= [ u{x)v{x)dx, ||M||' = («,n) (3.1) 

Jn 

for = (—L,L), respectively. Let Sn be the space of trigonometric polyno¬ 
mials of degree N/2 dehned as 

Sn = span{e*''"^/^| - N/2 < k < N/2 - 1} (3.2) 


where iV is a positive integer. : L^(fl) —)■ Sn is an orthogonal projection 
operator 

Ar/2-l 

Pnu{x)= (3.3) 

k=-N/2 

such that for any u G 

(m — P/vM, <p) = 0, V<p e Stv- (3.4) 
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The projection operator P^r commutes with derivative in the distributional 
sense: 

D^Pnu = PnD^u and = PnD]^u. (3.5) 

Here and P" stand for the nth-order classical partial derivative with 
respect to x and t, respectively. denotes the periodic Sobolev space 

equipped with the norm 

OO 

l|w|lH 5^ (1 + (3-6) 


where “fc = 777 / The Banach space Xg = C^{[0,T]; 

2T Jq 

is the space of all continuous functions in Pp(hl) whose distributional deriva¬ 
tive is also in Pp(r2), with norm IImIIx^ = max (||M(t)||g-|-||M 4 (t)||^). In order 

ts [ 07 ] 

to prove the convergence of semi-discrete scheme, we need following lemmas. 


Lemma 3.1. 

that 


IM, llJy For any real 0 < fi < s, there exists a constant C such 


\\u - PtvmIU < Vm e p;(p). (3.7) 

Lemma 3.2. fl^l Assume that f G C'^(M), u,v E fl L°°{Q) and 

k = [s] -|- 1, where s > 0. Then we have 


\\f{u)-f{v)\U<C{M)\\u-v 


(3.8) 


if ll'i^lloo < Al, ||n||oo < M, ||m||s < M and ||n||s < M, where C{M) is a 
constant dependent on M and s. 


Corollary 3.3. Assume that f G C'^(M) and u,v E P^(P) n L°°(f2) then 

\\f{u) - f{v)\\ 2 <C\\u-v \\2 (3.9) 

where the constant C depends on ||m||oo, IIi^IIoo and ||m|| 2 , lli^lb- 

3.2. Convergence of the semi-discrete scheme 

The semi-discrete Fourier Galerkin approximation to fll.3p - fll.4p is 

Utt = + VlU^xtt - mUxxxxtt + PNf{u^)xx, (3.10) 

u^{x,0) = PN(j){x), uf^{x,0) = Pp^'ifix) (3.11) 

where u^{x, t) E for 0 < f < T. We now state our main result. 
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Theorem 3.4. Let s >2 and u{x,t) be the solution of the 2L periodic initial 
value problem fll.3p - fll.4p satisfying u{x, t) E C^([0, T]; Hf{Ll)) for any T > 0 
and {x,t) be the solution of the semi-discrete scheme H3.10\) - l[3.11\) . There 
exists a constant C, independent of N, such that 

\\u-u^\\x,<C{T,r],,T]2)N^-^\\u\\x, (3.12) 

for the initial data E 

Proof. Using the triangle inequality, it is possible to write 

\\u — u^\\x2 < ~ Pnu\\x2 + \\Pnu — u^\\x2- (3.13) 

Using Lemma 13.11 we have the following estimates 

\\iu-Pxu)it)\\2<CN^-muit)\\s 


and 

\\{u-Pxu)t{t)\\2<CN^-^\utms 


for s > 2. Thus, the estimation of the hrst term at the right-hand side of the 
inequality fl3.13p becomes 


u — Pxu\\x2 ^ CN'^ ^|| m || x „- 


(3.14) 


Now, we need to estimate the second term \\Pxu — u^\\x 2 at the right-hand 
side of the inequality 03.131) . Subtracting the equation 03.1Up from 01.3p and 
taking the inner product with (p E Sx we have 

(3.15) 


Since 


(("U ((^ xxi ((^ PN'^^xxtti 9^) ((^ Pn’^^xxxxU') 9^) 0 

(3.16) 


for all G S'm, by 03.51) . the equation 03.15P becomes 

{{{Pnu - u^)tt - (Pnu - u^) XX - PliPxU - U^)xxU + V2{PnU - U^) xxxxtt 

- (/(w) -f{u^))xx}, P) =0 


(3.17) 
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for all G S'at. Setting tp = (P/vm — u^)t in ll3.17p . using the integration by 
parts and the spatial periodicity, a simple calculation shows that 

{{Pnu - u^)tt, {Pnu - u^)t) = ~ ( 3 - 18 ) 

{{Pnu-u^)^^, {PNU-u^)t) =-^^\\{PNU-U^)x{t)\\‘^, (3.19) 

{{Pnu-u^Uu, {Pnu - u%) =-~\\{Pnu - (3-20) 

{{Pnu - 

^ ^xxxxtt-) {P„u - u"),) = ~J(PnU - «''W(i)ir. ( 3 . 21 ) 

Substituting (I3.18p - (l3.2ip in (I3.17p . we have 

~ [||(Pa.« - «'"),(i)|p + ||(Pa.« - + vilKP^^u - -"(..(idt 

+ fj 2 ll(PNU - U^)a;xt{t)\\‘^] = {{f{u) - /(m^))xx, {PnU - U^)t) ■ 

(3.22) 


In the following, we will estimate the right-hand side of the above equation. 

Using the Cauchy-Schwarz inequality and the Corollary 13.31 we have 

{{f{u) - f{u^))^^, {PNU-U^)t) < \\{f{u) - f{u^))^^\\ \\{PNU-U^)t{t)\\ 

< l(ll/t“)-/(«~)ll^+ !!(/’»« 

(3.23) 


Substituting fl3.23p in fl3.22p we have 

+h2||(^ivn - u^)^^t{tW) 

<C{\\{u-u^){t)\\l + \\{PNU-u%{t)r) 

< C (||(« - PNumWl + ||(P^« - «^)(t)||^ + \\{Pnu - n^)*(t)|p) . 

(3.24) 


Adding the terms [Pnu — , {Pnu — u^)t) and {{Pnu — u^)xxi {Pnu — u^)xxt) 

to both sides of the equation fl3.24l) becomes 


1 

2 


min{l,771,772}^ \\{Pnu - u^){t)\\l + \\{Pnu 


u% 


<C{\\{u- PNumWl + \\{Pnu - u^mWl + \\{Pnu 


-u%{t)\\l). 
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Note that ||(P/vm — m^)(0)II 2 = 0 and IKP/v^ — M'^)t(0)||2 = 0. The Gronwall 
Lemma and Lemma 13.11 imply that 

\\{PNU-u^){t)\\l + \\{PNU-u^)t{t)\\l < f ||M(r) - PjvM(r)||2 

Jo 

< C{T,T]i,T]2) f ||M(r)||2dr 

Jo 

(3.25) 


for s > 2. Finally, we have 


\\PNU-u^\\x,<CiT,r],,r]2)N^-^ ||n|U.. 


(3.26) 


Using fl3.14p and (13.261) in fl3.13p . we complete the proof of Theorem 2.4. □ 

Note that the convergence theorem is proved for the semi-discrete Fourier 
Galerkin approximation. We point out that the Fourier collocation pseudo- 
spectral method is used in the following section as it is more practical due 
to the use of FFT. 


4. The fully-discrete scheme 

We solve the HBq equation by combining a Fourier pseudo-spectral method 
for the space component and a fourth-order Runge Kutta scheme (RK4) for 
time. If the spatial period [—L, L] is, for convenience, normalized to [0, 27r] 
using the transformation X = 7r(x + L)/L, the equation fl3.10p becomes 


u 


N 

tt 



N 

Uxx-Vi 



Uxxtt+V2 



4 

7 /^ - 

^XXXXtt — 



(4.1) 


The interval [0, 27r] is divided into N equal subintervals with grid spacing 
XX = 2n/N, where the integer N is even. The spatial grid points are given 
by Xj = 2'Kj/N, j = 0,1,2,N. The approximate solutions to u^{Xj, t) is 
denoted by Uj{t). The discrete Fourier transform of the sequence {Uj}, i.e. 


N-l 


1 N N 

Ut = T,\U,] = -Y,Ujm-i>{-ikX,), - <k<--l. (4,2) 


1=0 










gives the corresponding Fourier coefficients. Likewise, {t/j} can be recovered 
from the Fourier coefficients by the inversion formula for the discrete Fourier 
transform fl4.2|) . as follows: 

Uj = J^~^[Uk]= ^ t/fcexp(zfcXj), j = 0,1,2,1 . (4.3) 

2 

Here J- denotes the discrete Fourier transform and its inverse. These 
transforms are efficiently computed using a fast Fourier transform (FFT) 
algorithm. In this study, we use FFT routines in Matlab (i.e. fft and ifft). 

Applying the discrete Fourier transform to the equation fl4.ip we get the 
second order ordinary differential equation which can be written in the fol¬ 
lowing system 


{Uu)t = 

Vk 

(4.4) 

( 14 ), = 

1 + r]i{7rk/Ly + r]2{7Tk/Ly ^ ^ 

(4.5) 


In order to handle the nonlinear term we use a pseudo-spectral approxi¬ 
mation. That is, we use the formula J^k[{UjY] to compute the Fourier 
component of u^. We use the fourth-order Runge-Kutta method to solve the 
resulting ODE system (I4.4p - (l4.5p in time. Finally, we hnd the approximate 
solution by using the inverse Fourier transform (14.3p . 

5. Validation of the fully-discrete scheme 

The purpose of the present section is to verify numerically that (i) the 
proposed Fourier pseudo-spectral scheme is highly accurate, (ii) the scheme 
exhibits the fourth-order convergence in time and (Hi) the scheme has spec¬ 
tral accuracy in space. The Loo-error norm is dehned as 

Loo-error = max \ Ui — Ui \ (5.1) 

i 

where Ui denotes the exact solution at u{Xi,t). 

To validate whether the Fourier pseudo-spectral method exhibits the ex¬ 
pected convergence rates in time we perform some numerical experiments. 
In this section we use quadratic nonlinearity where f{u) = v?. We take hxed 
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number of spatial grid points N and various values for the number of tem¬ 
poral grid points M. For f]! = 'r ]2 = ^, the initial data corresponding to the 
solitary wave solution fl2.6p for quadratic nonlinearity become as follows: 

■u(a;,0) = Asech‘^{B x), (5.2) 

v{x,0) = 4:AB c sech"^ {B x) tanh {B x). (5.3) 

This initial data generate a solitary wave initially at Xq = 0 moving to the 
right with the amplitude A k, 0.39, speed c ~ 1.13 and B 0.14. The 
problem is solved for times up to T = 5. The space interval is chosen as 
— 100 < X < 100 to eliminate the error due to boundary effects. In these 
experiments, we take = 512 to ensure that the error due to the spatial dis¬ 
cretization is negligible. The convergence rates calculated from the Loo-errors 
at the terminating time T = 5 are shown in Table I. The computed conver¬ 
gence rates agree well with the fact that Fourier pseudo-spectral method 
exhibits the fourth-order convergence in time. 

TABLE I 


The convergence rate in time calculated from the Loo-errors in the case of single 
solitary wave {A ~ 0.39, N = 512). 


M 

Loo-error 

Order 

2 

8.662E-3 

- 

5 

2.530E-4 

3.8561 

10 

1.614E-5 

3.9704 

50 

2.623E-8 

3.9903 

100 

1.637E-9 

4.0021 


To validate whether the Fourier pseudo-spectral method exhibits the ex¬ 
pected convergence rate in space we perform some further numerical experi¬ 
ments for various values of N and a hxed value of M. In these experiments we 
take M = 1000 to minimize the temporal errors. We present the Loo-errors 
for the terminating time T = 5 together with the observed rates of conver¬ 
gence in Table II. These results show that the numerical solution obtained 
using the Fourier pseudo-spectral scheme converges rapidly to the accurate 
solution in space, which is an indicative of exponential convergence. 
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TABLE II 


The convergence rates in space calculated from the Loo-errors in the case of the 
single solitary wave {A ss 0.39, M = 1000). 


N 

Loo-error 

Order 

10 

0 .211E-1 

- 

50 

1.747E-3 

1.5480 

100 

4.431E-7 

11.9450 

150 

6.500E-10 

16.0916 

200 

3.884E-13 

25.8017 


6. Numerical Experiments 

In this section, we investigate the effect of extra dispersion term and 
the nonlinear term on the solntions of the HBq eqnation. For this aim, we 
will consider propagation of a single solitary wave, head-on collision of two 
solitary waves and blow-np solntion. 

6.1. Propagation of a Single Solitary Wave 

First, we study the effect of the nonlinear term on the single solitary 
wave solution. The problem is solved on the space interval —100 < x < 100 
for times up to T = 5 by using the initial data corresponding to the solitary 
wave solution fl2.6p . We show the variation of L^-eiiois with N for the 
HBq equation for various powers of nonlinearity, namely, f{u) = for p = 
2,3,4, 5 in Figure 1. The value of M is chosen to satisfy z/ = At/Ax = 
2.56 X 10“^. We observe that the Loo-errors decay as the number of grid 
points increases for various degrees of nonlinearity. Even in the case of the 
quintic nonlinearity, the Loo-errors are about 10“^^. We have not seen any 
numerical results for the HBq equation in the literature to compare with our 
results. This experiment shows that the proposed method provides highly 
accurate numerical results even for the higher-order nonlinearities. 

The question arises naturally how the solutions of the HBq equation be¬ 
have when the coefficient of extra dispersion term r ]2 0^. For this aim, 

we perform some numerical tests for various values of ri 2 and the hxed value 
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N 


Figure 1: Lqo- errors for the increasing values of N. 


rji = 1 for quadratic nonlinearity. We use the initial data 


u{x, 0) = y4sech^(- 


X) 


v{x, 0) = 2 ydsech^ 




( 6 . 1 ) 

( 6 . 2 ) 


where y4 = 0.9 and c = y —h 1. This initial data correspond to the initial 

prohle of the exact solitary wave solution for the IBq equation given in 0. 
In Figure 2, the solid line shows the exact solution of the improved Boussinesq 
equation which actually corresponds to the HBq equation with r ]2 = 0- The 
dashed lines show the numerical solution of the HBq equation for the values 
of r ]2 = 10, rj 2 = 5, rj 2 = 1 and r ]2 = 0.1 and rji = 1. The problem is solved 
on the space interval —100 < x < 100 for times up to T = 5 where N = 512 
and M = 5000. To observe the difference of two solutions more clear, this 
hgure is presented on the space interval —30 < x < 30. The numerical tests 
indicate that, as the parameter r ]2 tends to zero, the solutions of the HBq 
equation converge to the solitary wave solution of the IBq equation. 
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Figure 2: The numerical solution of the HBq equation for rj 2 = 10, r ]2 = 5, ri 2 = 1, 
r ]2 = 0.1 with r]i = 1 and the exact solitary wave solution of the IBq equation 


TABLE III 

The amplitude at the final time T = 5 with respect to changing r] 2 . 


h2 

Amplitude (p = 2) 

Amplitude (p = 3) 

0.1 

0.894 

0.943 

0.3 

0.885 

0.930 

0.5 

0.876 

0.919 

0.8 

0.870 

0.912 

1 

0.866 

0.908 


Table III shows the dependence of the amplitude to a changing ri 2 . In these 
experiments, we use the initial data (I6.ip - (l6.2p with N = 512 and M = 5000. 
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We note that initial amplitude is 0.9. From Table III, we can draw two 
observations: 

(i) The amplitude of the wave is decreasing with increasing values of 772 . 
This numerical result agrees well with the fact that the wave spreads out 
with increasing dispersive effects. 

(ii) The amplitude is increasing with increasing values of power of nonlinear¬ 
ity p. This numerical result is also agrees well with the fact that the wave 
steepens with increasing nonlinear effects. 

6.2. Head-on Collision of Two Solitary Waves 

In this section, we consider the HBq equation with quadratic nonlinearity 
and we study the head-on collision of two solitary waves with equal ampli¬ 
tudes. The initial conditions are given by 
2 

u{x, 0) = ^ As,ec\i^{B{x — Xg)), 

i=\ 

2 

v{x, 0) = 4 ^ ABciSech.^ [B{x — Xq)) tanh (B(x — Xq)) . 

i=l 

We consider two solitary waves, one initially located at Xq = —40 and moving 
to the right with amplitude A (ci > 0) and one initially located at Xq = 
40 and moving to the left with amplitude A (c 2 = —Ci). The problem is 
solved again on the interval —100 < x < 100 for times up to T = 72 
using the Fourier pseudo-spectral method. All experiments in this section 
are performed for = 512 and M = 7200. We hrst illustrate the head-on 
collision of two solitary waves with equal amplitudes A ^ 0.39 with rji = 
?72 = 1 and A ^ 1.08 with 771 = 772 = 2 in Figure 3. 

As the HBq equation cannot be solved by the inverse scattering method, 
the interaction of solitary waves are inelastic. Although secondary waves 
exist in all nonlinear interactions, they become more visible as we increase 
the amplitudes of the interacting waves. 

Since an analytical solution is not available for the collision of two soli¬ 
tary waves, we cannot present the Loo-errors for this experiment. But, as 
a numerical check of the proposed Fourier pseudo-spectral scheme, we have 
observed the evolution of the change in the conserved quantity Ji. The error 
I /i(t) — /i( 0 ) I is approximately of order 10 “^^ and 10 “^° for the amplitudes 
A ^ 0.39 and A ^ 1.08 up to the hnal time T = 72, respectively. This 
behavior provides a valuable check on the numerical results. 
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Figure 3: Surface plot of the head-on collision of two solitary waves with equal 
amplitudes. A « 0.39 (left panel) and A « 1.08 (right panel). 


6 .S. Blow-up Solution 

In the last section, we investigate the effect of extra dispersion term and 
the nonlinear term on the blow-up solutions of the HBq equation. The HBq 
equation is one of the class of nonlocal equation, we refer to Theorem 5.2 in 


10| as the blow-up criteria. Applied to the HBq equation, this theorem can 


be restated as: 


Theorem 6.1. Suppose cf) = = Tj, for some d>, T G and F{lp) E 

If there are some /r > 0 such that 

uf {u) < 2/iu^ -|- 2 (1 -|- 2/i) F (u) for all u G M, (6.3) 


and I 3 (0) < 0 , then the solution u of the Cauchy problem for HBq equation 
with the initial conditions m(x,0) = 0(a;), Ut{x,0) = fj{x) blows up in finite 
time. 


First, we consider the HBq equation with both quadratic nonlinearity f{u) = 
and cubic nonlinearity f{u) = —u^. We set the parameters r/i = 772 = 1 in 
fll.311 . For the quadratic nonlinearity, if we choose the initial data 

\_ jP 

— 1 e 3 ^ 'f)[x) = 


— l) e 2 (6.4) 


/2x^ 

= 4 — 
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Figure 4: The variation of ||?7||oo with time for increasingly refined discretization 
for quadratic nonlinearity (left panel) and cubic nonlinearity (right panel). 

as in [l^, and the condition (16.dh is satisfied for /i = Similarly, for the 
cnbic nonlinearity, the initial data are chosen as 

0(x) = 13 e il){x) = (l — e ^ (6.5) 

where the ineqnality fl6.3p is satisfied for /i = |. We note that /3(0) < 0 for 
both cases. For the experiments in this section, the problem is solved on the 
interval —10 < x < 10 np to hnal time T = 4 for qnadratic nonlinearity and 
T = 0.4 for cubic nonlinearity. 

In order to check whether the numerical blow-up results are affected by 
the time and space step sizes, we now make some numerical experiments 
with varying resolutions. We start with the spatial and temporal resolutions 
corresponding to = 64 and M = 500, respectively. Then we refine the 
mesh by increasing the number of both spatial and temporal grid points. We 
present the variation of the Loo-norm of the approximate solution obtained 
using the Fourier pseudo-spectral scheme for various mesh in Figure 4. To 
observe the blow-up profile more clear, figures are illustrated near the blow¬ 
up time. It is seen that curves converge to a limiting blow-up profile as the 
resolution is increased. For the number of spatial and temporal grid points 
more than N = 512 and M = 4000, the blow-up prohle is not distinguishable 
from each other. Therefore, we set = 512 and M = 4000 for the rest of 
the study. The numerical results strongly indicate that a blow-up is well 
underway by time t = 3.8 for quadratic nonlinearity and t = 0.35 for cubic 
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nonlinearity. In Figure 5, we also illustrate the blow-up profile near the blow¬ 
up time at t = 3.710, t = 3.711, t = 3.7115 and t = 3.712 for the quadratic 
nonlinearity as in [l^. We observe that the Loo-norm of the blow-up solution 
increases very rapidly while the blow-up prohle remains the same. Therefore, 
the blow-up appears to be a similarity type. 




X X 




X X 

Figure 5: The blow-up profile of the HBq equation near the blow-up time at t = 
3.710, t = 3.711, t = 3.7115, t = 3.712 for quadratic nonlinearity with rji = r ]2 = 1- 

To investigate how the blow-up time depends on the coefficient of extra 
dispersion term 772 with the hxed value pi = 1, we perform some numerical 
experiments. The variation of the blow-up time as 772 —t O’*" is presented in 
Figure 6. As it is seen from the hgure, less dispersion gives rise to earlier 
blow-up time. 

To investigate how the blow-up time depends on the power p of nonlinear 
term, we perform some numerical experiments. We use the initial data fl6.4l) 
for f{u) = uP, p = 2,4,6 and (16.51) for f{u) = —p = 3,5. We set the 
parameters 771 = 772 = 1 in (11.31) . We note that it is possible to hnd some 
appropriate values for p to satisfy the condition (16.31) for all these cases. 
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Figure 6: The variation of the blow-up time as rj 2 —)• O"*" for quadratic nonlinearity 
(left panel) and cubic nonlinearity (right panel). 




0.1 0.2 0.3 


t 


Figure 7: The blow-up time with respect to changing p. 


Moreover, 13 ( 0 ) also remains negative for increasing p with the initial data 
fl6.4p and fl6.5p . Figure 7 shows the blow-up time with respect to changing 
p. It can be observed that the blow-up time is decreasing with increasing 
powers of nonlinearity. 
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